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ABSTRACT 


Some  simple  transformation  properties  of  the  diffusion 
equation  that  has  been  classically  used  to  describe  turbulent 
transport  in  the  atmosphere  sre  briefly  examined.  The  essen¬ 
tial  ten8orial  nature  of  the  eddy  diffusivity  1b  emphasized, 
and  it  is  concluded  that  the  standard  form  frequently  adopted 
for  the  diffusion  equation  in  the  meteorological  literature 
cannot  be  generally  valid. 
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I.  INTRODUCTION 


We  consider  the  equation  that  is  frequently  quoted  in  the  meteorologi¬ 
cal  literature  as  describing  the  turbulent  diffusion  of  a  conservative 
scalar  property  in  three  dimensions.  For  present  purposes  the  density  of 
the  fluid  medium  will  be  supposed  to  be  constant,  and  the  standard  general 
form  of  the  equation  taken  to  be1- 
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relative  to  a  fixed  system  of  rectangular  Cartesian  coordinates  Oxyz. 

Here  C  ?  C  (x,y,z,t)  is  a  suitably  defined  mean  concentration;  u,  v,  w 
are  the  mean  velocity  components  of  the  fluid  motion;  R^,  K2,  K-j  are  the 
so-called  eddy  diffusion  coefficients,  which  are  assumed  to  be  functions 
of  position.  Unfortunately  the  precise  nature  of  the  quantities  Kj,  K2, 
and  (e.g„,  whether  they  are.  scalar  functions  of  position  or  components 
of  i  vector  ete.l  is  rarely  discussed  in  the  meteorological  con¬ 

texts.  Since  this  question  is  vital  to  the  interpretation  or  rnu  equation 
it  is  briefly  examined  in  this  manuscript.  It  Is  shown  that,  even  under 
the  customary  assumptions.  Equation  (1)  cannot  be  a  generally  valid  form 
of  diffusion  equation. 


II .  A  GENERAL  TRAM3FOBMATIOM 


More  fundamentally,  the  right  side  of  Equation  (1)  is  the  negative  of 
div  F  s  SEx/ax  +  ^Fy/dy  +  &F*/az  where  £  -  (Fx,  Fy,  Fz)  is  the  turbulent 
flux  vector  for  the  scalar  property  considered.  Equation  (1)  then  results 
from  Che  hypothesis  that  the  components  of  F  relative  to  the  axes  Oxyz  are 
related  In  a  simple  homogeneous  linear  fashion  to  the  respective  components 

of  the  concentration  gradient  vector  ^^/br  3  QC/bx,  ^C/dy, 

Let  the  Greek  symbols  a,  P,  Y,  etc.  stand  for  any  of  x,  y,  or  z,  and  let 
us  define  for  the  chosen  axes  Oxyz  an  ordered  array  of  nine  quantities 
Kop  by  the  matrix  equation 
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For  convenience  we  shall  also  use  the  Einstein  summation  convention, 
according  to  which,  in  any  term  containing  a  repeated  Greek  suffix,  it 
is  understood  that  the  suffix  is  given  all  possible  values  x,  y,  z  and 
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the  terms  so  obtained  are  rhen  added.  Then  at  any  chosen  point  of  the 
medium  Equation  (1)  is  exactly  equivalent  to  the  following  relation 
between  the  components  of  F  and  those  of 
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We  now  consider  a  second  arbitrary  set  of  rectangular  axes  Ox'y'z',  and 
let  ta^  i  denote  the  cosine  of  the  angle  between  the  axes  Oci  and  Op',  so 
that  t^Qi  =  t_  t  *  Then  the  components  of  the  vector  F  relative  to  the  new 
axes  will  be  giver  by 


V  *=  Fa  taY' 

dC  / 

Similarly,  the  components  of  the  gradient  vector  jor  for  the  axes 
Oxyz  will  be  related  to  those  for  the  axes  Ox'y'z'  by- 
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From  Equations  (3),  (4),  and  (5)  it  follows  that 


V  =  -  KoP 


dC 


3rp  taY'  "  "  K«.g 


3c  I 


'P  a.Y' 


dC  / 

"  KaP  /^rj  »  taYl  tp6 ' 
rela 

7or  in  the  new  axes  Ox'y'z'.  It  may 

3C  / 

fY'  =  —  ’S'S  '  p r& ,  (7) 


Equation  (6)  provides  the  rule  for  relating  the  components  of  the  flux 
vector  F  and  the  gradient  vector  —  <"  n».i  avaa  iVn'.1  tc  mo. 
equivalently  be  written 


where 


N'S  '  =  KoP  taY'  UP6 ' 


(8> 


Equation  (8),  however,  defines  the  known  transformation  law  for  a  second- 
order  Cartesian  tensor.  It  thus  establishes  the  fact  that  the  set  of 
quantities  Kfjp  in  Eqviation  (2)  constitute  the  components  of  a  second- 
order  l^ensor.  This  is  a  direct  consequence  of  the  vectorial  nature  of 
F  and  °c/dr.  The  tensor  Ka^  is,  of  course,  in  general  a  function  of  posi¬ 
tion,  since  the  relation  be'tween  F  and  d(ygr  will  be  defined  at  every 
point  of  the  medium.  We  see  from  Equation-^)  that  in  the  arbitrary 
axes  Ox'y'z'  the  components  of  F  are  general  homogeneous  linear  functions 
of  the  components  of  ^/^r-  This  relati on ship  la  exactly  analogous  to 
that  encountered  in  consTder ing  molecular  heat  conduction  in  an  anisotropic 
medium.  Whether  such  a  linearity  hypothesis  is  actually  justified  for  the 
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turbulent  flux  of  a  scalar  property  can,  of  course,  only  be  decided  on 
the  basis  of  experiment.  For  any  physical  meaning  to  be  attached  to  the 
tensor  Ky'g'  (or  KofJ),  it  is  evident  that  it  must  be  a  property  of  the 
fluid  motion,  i.e.,  independent  of  the  distribution  of  the  diffusing 
scalar  property.  We  may  call  the  eddy  diffusivity  tensor, 

'ft  may  oe  noted  tWitt  the  requirement  for  such  a  tensor  quantity  is 
actually  suggested  by  an  obvious  extension  of  the  Prandtl  mixing-length 
concept  to  three  dimensions.  Thus,  if  we  suppose  that  the  turbulent 
fluctuation  of  concentration  C'  from  the  mean  value  C  is  given  by 


so  that  _1  =  (lx,  ly,  iz)  is  the  mixing-length  vector,  then  if  Va'  denotes 
a  component  of  the  turbulent  velocity  fluctuation,  and  an  overbar  denotes 
a  mean  va  lue , 

f  =  v'  cT  =  -  v1  i„ 

o  a  a  p  /or  g 

=  ~K  /  where  K  =  v!,  lft 

ap  I  dTp  aP  a  P 

and  K  ,,  is  clearly  a  second-order  tensor. 

up 

Relative  to  the  axes  Ox'y'as1  the  right  side  of  Equation  (1),  which 
represents  — div  F,  is  replaced  by  the  nine-term  sum  represented  under 
Che  summation  convention  by 


III.  SOME  SPECIAL  FORMS 


The  special  diagonal  form  of  the  tensor  in  Equation  (2)  relative 
to  the  axes  Oxyz  implies  that  at  the  chosen  point  of  the  medium  the  dif- 
fusivity  tensor  is  symmetric  (since  symmetry  is  a  property  preserved  under 
transformation  of  axes,  and  Kqo  is  clearly  symmetric)  and  further  that 
Oxyz  are  its  principal  axes.  For  Equation  (1)  to  be  valid  at  all  points 
of  the  medium  it  would  also  be  necessary  for  the  tensor  to  have  Oxyz  as 
principal  axes  at  all  points  of  space.  It  therefore  follows  that  Equation 
(1)  cannot  be  adopted  as  a  possible  general  form  of  diffusion  equation 
because  the  axes  Oxyz  cannot  be  chosen  arbitrarily  but  must  be  a  preferred 
set  with  the  above  properties.  A  special  case  would  arise  if  the  (symmetric) 
diffusivitv  tensor  were  isotropic  at  all  points  of  the  medium,  i.e.,  if 
Kf  =  K2  =  K3  =  say  k  (a  scalar  function  of  position)  so  that 

1  0  0  \ 
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In  this  case  its  components  would  be  the  same  for  all  rectangular  axes 
and  all  axes  would  be  principal  axes.  Then  from  Equation  (3),  the  flux 
vector  would  be 


F 


and  the  corresponding  diffusion  equation 
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For  this  special  isotropic  form  of  diffusion  equation  there  are  no  pre¬ 
ferred  axes. 

Having  recognized  the  above  fundamental  limitations  on  the  use  of 
Equation  (1)  we  may  seek  for  a  meteorological  situation  where  it  might 
appear  plausible  to  postulate  the  existence  of  a  preferred  set  of  axes, 
independent  of  spatial  location.  Such  a  situation  could  occur  in  the 
surface  layers  of  the  atmosphere,  where  the  mean  wind  vector  can  be 
regarded  approximately  as  everywhere  parallel  to  a  given  vertical  plane. 
This  uniquely  defined  vertical  plane,  together  with  the  vertical  direction 
(which  must  be  fundamental  to  the  characteristics  of  turbulence  generated 
by  the  horizontal  ground  surface)  will  suffice  to  define  a  preferred  set 
of  axes,  e.g.,  Ox  and  Oz  in  the  plane  and,  respectively,  horizontal,  and 
vertical;  and  Oy  perpendicular  to  the  plane,  i.e.,  perpendicular  to  the 
mean  velocity  vector.  Relative  to  these  axes  the  mean  velocity  component 
v  of  Equation  (1)  is  v  ^  0.  If  it  is  postulated  that  this  preferred  set 


of  axes  is  a  principal  set  for  a  symmetric  diffusivity  tensor  at  every 
point  of  the  medium,  then  relative  to  this  special  set  of  axes  the  dif¬ 
fusion  equation  would  have  the  form 
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where  Kj,  K2/  and  K3  are  soalar  functions  of  position  that  define  the 
principal  components  of  the  diffusivity  tensor.  For  the  situation 
visualized  above,  it  is  also  normally  assumed  that  the  vertical  compon¬ 
ent  of  mean  velocity  w  =  0.  Of  course,  whether  auch  a  set  of  stringent 
postulates  does  correspond  to  the  real  atmospheric  diffusion  process  can 
only  be  justified  on  the  basis  of  experiment. 


It  should  be  emphasized  that  the  possible  existence  of  a  preferred 
set  of  axes  in  the  above  example  depends  on  the  assumption  that  the  mean 
velocity  vector  is  everywhere  parallel  to  a  given  vertical  plane,  so  that 
a  preferred  set  of  horizontal  axes  can  be  defined  (since  we  have  already 
assumed  that  the  vertical  direction  is  a  preferred  one).  If  this  assump¬ 
tion  were  not  made,  for  example,  if  we  were  to  assume13  only  that  the  mean 
velocity  was  parallel  to  the  ground  surface  with  general  components 
u  f  0,  v  f  0,  w  “  0,  then  no  preferred  horizontal  axes  Ox,  Oy  could  be 
defined,  and  hence  a  diffusion  equation  of  the  form  (K^  f  K2) 
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would  not  be  va  lid,  This  form  of  equation  could  only  be  valid  in  the 
special  case  Kj  =  K2,  i.e.,  if  the  tensor  quadric  had  rotational  symmetry 
about  the  z  axis. 

Finally,  it  should  be  noted  that,  although  the  above  discussion  has 
been  restricted  to  rectangular  axes  for  reasons  of  simplicity,  the  ques¬ 
tions  raised  will  be  quite  fundamental  when  considering  possible  forms 
for  the  equation  of  turbulent  diffusion  in  other  coordinate  systems,  e.g., 
any  system  of  orthogonal  curvilinear  coordinates.  For  such  cases  it  will 
be  necessary  to  apply  Che  methods  of  general  tensor  analysis. 
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